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ABSTRACT 

Global time transfer, or synchronization, between a user 
clock and USNO UTC time can be performed using the Global 
Positioning System (GPS), and commercially available time 
transfer receivers. This paper presents the test results 
of time transfer using the GPS system and a Stanford Tele- 
communications, Inc. (STI) Time Transfer System (TTS) 

Model 502, Tests at the GPS Master Control Site (MCS) in 
Vandenburg, California and at the United States Naval 
Observatory (USNO) in Washington, D.C. are described. An 
overview of GPS, and the STI TTS 502 is presented. A 
discussion of the time transfer process and test concepts 
is included. 


INTRODUCTION 

This paper describes test results using the Stanford Telecommunications Inc. 
(STI) Time Transfer System (TTS) Model 502 to perform time transfer between 
a user clock and Universal Time Coordinated (UTC) available from the United 
States Naval Observatory (USNO). The time transfer procedure makes use of 
the STI TTS 502 receiver and the Global Positioning System (GPS) to synchro- 
nize a user's clock to UTC in a two step procedure:, 

1) time transfer between the user's clock and the GPS 
Master Clock (located at Vandenburg, California) 
using observations of GPS satellites, 

and 

2) transformation between the GPS Master Clock and 
UTC time based on USNO synchronization data. 

A brief description of the GPS system is provided along with a discussion 
of the time transfer procedure used at the USNO to synchronize UTC and GPS 
MCS time. More detail is available in the literature (References (1) and 
(2), respectively) . A description of the time transfer procedure using the 
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STI TTS 502 receiver, is also presented. 


Test results are described which indicate the absolute accuracy of time 
transfer using the STI TTS 502. These results are based on tests performed 
at the GPS MCS. Further test results at the USNO establish the baseline 
accuracy with which time transfer can be performed between a user clock 
and UTC. 


BACKGROUND 

The GPS system and commercially available receiver systems, provide the 
capability for a user to determine his global position, velocity and time 
with a high degree of accuracy. As originally planned, the GPS was to 
consist of a space segment of twenty- four satellites and a ground segment 
of a Master Control Site (MCS) and five or more Monitor Sites (MS), one of 
which was to be located at the USNO. The monitor sites function is to 
receive transmissions from each GPS satellite, referred to a local clock, 
and to retransmit this information to the MCS over secure data communications 
links. The MCS correlates this information with other data, performs the 
necessary calculations to determine current satellite performance parameters 
and uplinks this information to each GPS satellite on a daily or, as required, 
basis. The uplinked data provides current information on clock performance 
and navigation data. The satellites were intended to be equally distributed 
in three orbital planes inclined to the Earth's equatorial plane by 63 
degrees with the ascending nodes of the orbital planes differing by 120 
degree intervals. 

Funding cutbacks have reduced the planned number of satellites to eighteen. 
This cutback will affect position/velocity determination, but have little 
effect on time transfer, since only one satellite needs to be visible to a 
user to perform time synchronization. 


GPS SATELLITE SIGNAL 

Reference (3) describes the GPS data transmission process and structure. 
Briefly, data are transmitted from the GPS sateTlites on two carrier fre- 
quencies, a primary, Lj, at 1575.42 MHz and a secondary, L 2 at 1227.6 MHz. 

The Li frequency is simultaneously modulated by a precisiori (P) code and a 
coarse/acquisition (C/A) code. The L? frequency is modulated by either a 
P or C/A code. Every six seconds a 50 bps message stream is transmitted 
with a total frame size of 1500 bits, common to both frequency bands and 
to both P and C/A codes. Each frame is divided into five 300 bit subframes 
which are further subdivided into ten 30-bit words. Telemetry and code 
handover information is contained in the first two words of each subframe. 

Key to the time transfer process is data contained in the first three sub- 
frames. The last eight words of Subframe 1 contain clock corrections con- 
sisting of second order polynomial coefficients which represent the variation 
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of the GPS satellite onboard clock with respect to the GPS master clock. 

Also included is an age of data word, and ionospheric delay model parameters. 
The GPS satellites ephemeris and age of data words are contained in Subframes 
2 and 3. Subframe 4 provides an alphanumeric message for GPS users while 
Subframe 5 contains an Almanac (reduced accuracy version of information in 
Subframes 2 and 3) for each of the other GPS satellites in the constellation. 
Subframe 5 contains data on a single GPS satellite so that a sequence of 
frames is required to specify the Almanac for the complete GPS constellation. 

The time transfer process using GPS requires that a user be capable of 
receiving, demodulating and decoding the transmitted data stream from a 
visible GPS satellite. The observed GPS satellite transmits the onboard 
satellite clock time (and corrections to the GPS master clock), which allows 
the user to compare the "observed" GPS time with time from his own clock. 
Corrections for propogation delays are required including free space 
instantaneous distance from GPS satellite to user, tropospheric/ionospheric 
delays, and equipment delays. Applying these corrections allows a measure- 
ment of the difference between the GPS master clock and the user's clock. 

The ST I TTS Model 502 performs such a measurement for a user at a known 
location. 


STI TIME TRANSFER SYSTEM 


The STI Time Transfer System is shown schematically in Figure (1) and 
described in detail in Reference (4), The system consists of a GPS antenna, 
a preamplifier, time transfer receiver/processor and alphanumeric terminal 
for control of the receiver/processor. An external user clock (e.g.. Cesium 
Beam Standard) is connected to the TTS 502 including a 5 MHz and 1 pps input. 
Firmware in the Motorola microprocessor provides several functions, including 
the following. 

1) Scheduling the receiver to observe user selected GPS 
satellites. An automatic update option allows a user 
to specify an initial observation schedule. The 
microprocessor then uses the GPS ephemeris data to 
compute look angles accounting for orbit precession. 

2) Perform raw time transfer observations every six 
seconds as scheduled, including; 

a) derive satellite transmission time, 

b) compute satellite clock error with respect 
to GPS system time, 

c) compute GPS satellite position at the GPS 
system time, 

d) estimate propogation delays (ionospheric, 
tropospheric, receiver constant bias), 

e) compute satellite to station range, 

f) determine the user clock error (difference 
between instantaneous user clock time and 
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corrected GPS time). 

3) Perform data smoothing for a single pass of raw 
time transfer measurements, 

4) Data display and formatting. 

5) Operator interface (via the televideo 920C terminal) 
for data base update (user location, receiver bias) 
and mode control . 

The TTS 502 thus allows a user clock to be compared with the GPS system 
time. The differences are output to the televideo 920C CRT screen and 
(optionally) to external peripherals via an RS-232 output channel. Cor- 
rection to UTC time may be performed by inputting linear slope and offset 
coefficients which represent the variation between GPS system time and UTC 
time. 


TIME TRANSFER CONCEPT 

The concept of time transfer using GPS satellites is illustrated in Figure 2. 
The GPS satellite system will ultimately consist of a space segment of 
eighteen satellites distributed in three orbital planes inclined to the 
Earth's equatorial plane, Each satellite contains onboard an atomic clock 
(Rubidium or Cesium). A GPS ground segment includes a Master Control Site 
(MCS) and five or more Monitor Sites (MS). The monitor sites receive trans- 
missions from each of the GPS satellites, referred to a local clock, and 
retransmit this information to the MCS over secure data communication links. 
The MCS correlates the data and performs calculations to determine current 
satellite performance parameters, including clock performance, and uploads 
data to each GPS satellite on a daily or as required basis. This process 
is shown schematically in Figure 2 as an upload/ download from the MCS to 
a GPS satellite on the left hand side of Figure 2. 

Each GPS satellite continually transmits a message which repeats every six 
seconds. The message is available to users with special GPS receivers and 
contains a time signal, correction parameters between the particular GPS 
satellite clock and the GPS master clock, propogation correction parameters 
for ionospheric delay, the satellites ephemeris and clock age data. Other 
data in the message concerns the ephemerides of all the other GPS satellites 
in the constellation. 

Daily, for each GPS satellite the U.S. Naval Observatory 's (USNO) prototype 
TTS, performs a comparison between the individual GPS. -satellite clocks and 
the USNO master clock (hereafter referred to as Universal Time Coordinated 
(UTC)). In addition, the USNO compares UTC and the GPS time system. This 
data (viz., the correction between UTC and GPS, and the correction between 
UTC and the GPS satellite clocks) is published by the USNO on a daily basis 
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via the USNO Time Service Series 4. This process is depiGted in Figure 2 
in the center of the figure. 

As the GPS satellite becomes visible to the user site» the on-site TTS 
receives the GPS signal and compares the GPS time (via the onboard GPS 
satellite clock and GPS correction parameters ) to the user's clocks. The 
TTS equipment and software perform a running comparison between the site 
clock and GPS time. The data from the USNO is then applied to correct the 
UTC time. 


TEST PROCEDURES 

During June 17 to 19 of 1982, tests of the STI TTS Model 502 were performed 
at the GPS MCS in Vandenburg, California and during August and September at 
the USNO in Washington, D.C. These tests established measurements of the 
fundamental accuracy of time transfer using the STI TTS 502 and the accuracy 
of the process of time transfer from a user clock to UTC time. 


GPS MCS TESTS 


In order to establish an accuracy baseline of the Time Transfer Unit, an 
accuracy test was performed at the GPS MCS. The test is shown schematically 
in Figure 3. Two TTS-502 receivers (serial #002 and 003) were set up in the 
GPS Master Site Integration Facility (MSIF) located at Vandenburg AFB. 

A USNO portable cesium clock (PC 837) was used as an intermediate standard 
for input to each TTS receiver. Time transfer observations were made via 
the time transfer TTS-502 receivers and the GPS system. These observations 
provided measurements of the offset between PC837 and GPS system time, as 
shown schematically as the output of a TTS receiver at the bottom of Figure 
3. The observations were scheduled in accordance with the visibility of the 
GPS satellites at the MCS and to obtain a variety of time transfer measure- 
ments under all conditions (e.g. low and high elevation angles, day and 
night observations, etc.). The TTS-502 receivers were initialized to collect 
10 minutes of data per visible GPS satellite. A round robin scheme provided 
for ten minute data collection on each GPS satellite (SV numbers 4, 5, 6, 8 
and 9) in sequence. The TTS-502 smoothed the raw time transfer data, which 
were observed every six seconds throughout the ten minute data collection 
span, using a running 20 point linear least squares fit. A total of 149, 
ten minute data sets were observed during June 17 through 19. 

During the same period, direct measurements of the offset between PC837 and 
GPS system time were obtained using the USNO 5328 counter, as shown in 
Figure 3. The counter measured, hourly, the offset between PC 837 and the 
Vandenburg Monitor Station (VMS) clock. Cable delays were measured by a 
5370 Time Interval Counter to an accuracy of about 4 nsec. The difference 
between the VMS clock and GPS system time is monitored by GPS MCS personnel, 
as described in reference (5). A transformation was derived between PC837 
and GPS system time using the USNO counter observations and the VMS-GPS 
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offset. The relation is: 


^^’*'^PC837-GPS (0.041147 nsec/min)(T - T^pQ^^) - 44531.7 nsec 

where TppQQj^ = day 170, 4^, 15"’, 0^ and T is the observation time. An 

uncertainty of + 12.8 nsec is associated with the transformation, due to the 
errors associated with measuring the cable delays, noise in the USNO 5328 
counter, and uncertainty in the GPS MCS determination of the relation between 
VMS time and GPS system time. The absolute error in the TTS-502 time 
transfer observations can be determined by an expression of the form, 

(TTS-502 Error) = (^t)pc837-GPS ' ^^^^TTS-502 observations 

where 

(At)yTs _502 observations observations indicated as the output 

of the TTS-502 receivers {shown schematically at the bottom of Figure 3). 


USNO VERIFICATION TESTS 

Verification tests of the TTS-502 receiver, serial number 003, were performed 
at the USNO on August 19 through September 28 of 1982. The tests provided 
a measurement of the accuracy with which time transfer can be performed 
between UTC and GPS using the TTS-502 and the GPS satenite system. 


The test compared the prototype TTU currently at the USNO with the TTS-502. 
In particular, the prototype USNO TTU is currently being used to perform 
time transfer between the GPS master clock to UTC time, and the individual 
GPS on-board satellite clocks to UTC time (Reference (3)).. This is done 
at the USNO by connecting a physical UTC time signal into the prototype TTU 
and allowing the TTU to perform the time synchronization measurement. The 
results of these measurements are published daily by the USNO in a Time 
Series Report #4. The TTS-502 was simply set up beside the prototype TTU 
at the USNO. A UTC time signal was input to the TTS-502 as a system clock. 
Then, a comparison was made between the GPS time synchronization performed 
by the prototype TTU and that performed by the TTS-502. This test s'imply 
compares the operation of the prototype and new TTU units. Excluding 
observational noise which results directly from each TTU (viz. receiver 
noise, etc.), the time transfer measurements should be nearly identical, 
considering that the observational conditions are nearly identical (i.e. , 
same GPS satellite, same elevation angle, same ionospheric delays, etc.) 

A total of 154 passes of data were obtained for four GPS satellites (SV# 5, 
6, 8 and 9). Each pass of data was smoothed by a linear least squares fit 
to result in a single smoothed time transfer observation. The analysis of 
these 154 smoothed observations is described in the next section. 
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TEST RESULTS 

The results of both the GPS MCS tests and the USNO tests Indicate that the 
TTS-502 receivers provide the capability to perform time transfer measure- 
ments to an absolute accuracy which is well within the 100 nanosecond budget 
for the receivers. 


GPS MCS TEST RESULTS 

Figure 4 shows a histogram of the data collected at the GPS MCS. The 
histogram includes all data collected during June 17 to 19, for receiver 
serial numbers 002 and 003, excluding GPS satellite number 4. Receiver 
#002 shows an absolute bias of -56.0 nsec with a standard deviation of 
9.8 nsec. Receiver #003 exhibits a bias of -51.0 nsec with a standard 
deviation of 15.0 nsec. As shown on Figure 4, the one sigma uncertainty 
of the zero point calibration is +13 nsec. 

Subsequent to the GPS MCS test, ST I discovered a calibration error in the 
TTS-502 due to the algorithm used to correct for quantization effects . The 
resulting recal ibration reduced the TTS-502 bias error by approximately 25 
nanoseconds. Note that the STI recalibration was nbt fortuitous based on 
knowledge of the GPS MCS tests . STI personnel did not have the results of 
the GPS MCS test when the recalibration was performed. Thus , the absolute 
bias in the STI TTS-502 receivers is in the range of -26 to -56 nanoseconds 
with observational noise of 9 to 15 nanoseconds . It should be noted that 
the MCS GPS tests were conducted in an uncontrol led environment with 
variations in temperature, atmospheric pressure and humidity. Moreover, the 
data used to obtain the noise and bias figures, included observations of all 
GPS satellites (excluding number 4), day and night observations, and high 
and low elevation angles. Control of the ambient environment and careful 
selection of observations should decrease the observation noise. 


USNO TEST RESULTS 

Receiver #003 was moved to the USNO for further testing duri ng August 19 
through September 28. ,A total of 154 passes of time transfer data were 
collected in accordance with the previously described USNO test procedures. 
The data consisted of time transfer measurements between UTC and GPS system 
time using both the TTS-502 receiver and the USNO prototype receiver. Three 
basic analysis were performed including; (1) a comparison of the TTS-502 
and USNO prototype receivers, (2) an analysis of the abi 1 i ty of a linear 
model to represent the TTS-502 time transfer data, and (3) an analysis of 
the accuracy with which the UTC-6PS offset can be predicted using a linear 
model . Each of these is described below. 

First, a comparison was made between the TTS-502 time transfer data and that 
observed by the USNO prototype. Based on the results of the GPS MCS tests 
it was expected that the comparison of the TTS-502 and USNO prototype would 
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yield an offset within a factor of two of the expected error of the TTS- 
502. 

A linear regression was perfonned between the TTS-502 data and the prototype 
data. The resulting linear relation is, 

(TTS-502 Offset) = ( 1.000) (Prototype Offset) + 170. nsec. 

Figure 5 shows a histogram of the residuals about the linear expression. 
Clearly, a constant bias of 170 nanoseconds exists between the TTS-502 and 
the prototype receivers. This is significantly larger than the measured 
(-51 nsec) bias of the TTS-502 receiver. The standard deviation of the 
residuals about the regression line was 8.1 nsec. 

During the USNO testing, another time transfer receiver from the Naval 
Research Laboratory (NRL) was brought to the USNO by Mr. Jay Oaks. Initial 
results indicate that the TTS-502 compares well with the NRL receiver. 
Further tests are required to resolve the offset issue. STI plans to 
recalibrate the prototype receiver during December 1982. 

The second analysis performed on the USNO data was to model the observed 
TTS-502 data by a linear model. The data were fit by linear regression over 
several data collection spans. Expressions were obtained of the form, 

(GPS - UTC Offset) = A(0bservation Time - Epoch) + B 

where the observation time is the starting time of the tracking of a GPS 
satellite in Julian days (format DDDD.DD). The epoch was Chosen as Modified 
Julian day number 45198.0. The quantity A is the slope of the variation of 
the GPS to UTC time scales (microseconds per day), and B is the offset 
(microseconds) between GPS and UTC time at epoch. 

Linear regressions were performed using all satellite data for several 
data collection spans ranging from a few days to the entire data collection 
span of 40 days. For each regression, the standard deviation of the 
residuals about the regression equation was computed. The results are 
shown in Figure 6 which plots the standard deviation of the residuals as a 
function of the data collection span. It might be expected that the 
standard deviation would be small for short data collection intervals, 
indicating a "local fit" phenomena, and subsequently increasing and leveling 
off for the data collection span for which the linear model is valid. At 
some point, the standard deviation of the residuals would increase when a 
data span is reached for which the linear model becomes invalid. Such a 
phenomena is shown in Figure 6. The linear model provides a good fit for 
a data collection span of less than 35 days. Beyond that interval, the 
linear model seems to break down. Because of the abrupt change of the 
standard deviation of the residuals, it might be suspected that some change 
in the GPS time scale was effected at Modified Julian day number 45233. 
Similar changes in the GPS system time scale have occured in the past. 
However, over a collection interval of up to 35 days, a linear relation 
provides a good model to represent the time transfer data. 
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The final analysis of the USNO data was an investigation of the ability of a 
linear model to predict the offset between GPS and UTC time. Specifically, 
the TTS-502 data was fit by a linear regression over a 26 day subinterval 
of the 40 day span of collected time transfer data. An equation of the form, 

(GPS - UTC Offset) = A (Observation Time - 45198.0) + B 

was obtained by least squares fit of the TTS-502 data for Modified Julian 
days 45198 to 45224. Five separate regression equations were obtained, one 
for each GPS satellite (numbers 5, 6, 8 and 9) and one regression which 
included all four satellites. Figure 7 lists the regression equation co- 
efficients, the number of observations included in the fit, and the standard 
deviation of the residuals about the regression equation. 

Each of these equations, or linear models for a 26 day span of time transfer 
observations, was used to predict the GPS-UTC offset for times beyond the 
end of the fitting interval (i.e., beyond Modified Julian day 45224). These 
predictions were compared against the actual TTS-502 observational data 
beyond day 45224. The comparison. 

Prediction Error = (Predicted Offset) - (TTS-502 Observed Offset) 

was made to determine how the linear fit degrades beyond the end of the fit- 
ting interval. Figure 8 graphs the prediction error versus prediction 
interval (observation time -45224) for all four GPS satellites. The error 
degrades quadratically with prediction interval. However, a linear model 
yields a prediction error of less than 100 nanoseconds for a prediction 
interval of four to five days, using a 26 day calibration span. 


GPS SV 
NUMBER 

NUMBER OF 
OBSERVATIONS 

REGRESSION EQUATION COEFFICIENTS 

STANDARD 

DEVIATION 

RESIDUALS 

(USEC) 

SLOPE (A) pSEC 

OFFSET (B) pSEC 

5 

17 

0.110 

-39.5 

0.013 

6 

27 

0.109 

-39.1 

0.017 

8 

27 

0.111 

-39.2 

0.041 

9 

19 

0.110 

-39.1 

0.019 

ALL 

90 

0.110 

-39.1 

0.027 


FIG. 7 REGRESSION EQUATIONS FOR A 26 DAY OBSERVATION SPAN 
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ABSOLUTE VALUE OF THE PREDICTION ERROR (NSEC) 
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QUESTIONS AND ANSWERS 


MR. D. W. ALLAN, National Bureau of Standards 

A question regarding the bias and the distribution of that bias: The 9.8 
and 15 nanosecond widths, were those calculated from the 10-minute averages 
for all of the data taken for a receiver? I'm not sure how that 9.8 nano- 
second width was calculated. 

MR. D. L. HALL, HRB-Singer, Inc. 

Yes. Basically, we used each data point as a smooth data point that con- 
tained for, say* a 10-minute span, that 9.8 

MR. ALLAN: 

So each data point was a l0-m1nute average? 

MR. HALL; 

That's right. 

MR. ALLAN: 

And that's the distribution of those 10-minute averages? 

MR. HALL: 

That's correct. 

MR. ALLAN: 

Okay. I think in principle, from the work we have done at the Bureau, we 
have seen distributions between receivers looking at the same satellite 
that are roughly on an order of magnitude better than that, and I wonder 
if you have any feeling for the source of that distribution? 

MR. HALL: 

No. I did want to say, though, that those numbers represented observations 
that were taken both at low elevations, high elevations, all times of the 
night and day, and so it is not the best case. When we looked at using, 
say, a single satellite, using night observations, high elevations, we 
could certainly reduce those noise figures by certainly 4 or 5 nanoseconds 
but not by an order of magnitude. 

MR. ALLAN; 

That certainly explains a lot of It. 
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